Effective coupling between two Brownian particles 
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We use the system-plus-reservoir approach to study the dynamics of a system composed of two 
independent Brownian particles. We present an extension of the well-known model of a bath of 
oscillators which is capable of inducing an effective coupling between the two particles depending 
on the choice made for the spectral function of the bath oscillators. The coupling is non-linear 
in the variables of interest and an exponential dependence on these variables is imposed in order 
to guarantee the translational invariance of the model if the two particles are not subject to any 
external potential. The effective equations of motion for the particles are obtained by the Laplace 
transform method and besides recovering all the local dynamical properties for each particle we end 
up with an effective interaction potential between them. We explicitly analyze one of its possible 
forms. 

PACS numbers: 05.40.Jc 



Usually the Brownian motion of a given dynamical 
variable is modelled by considering the system it de- 
scribes - the system of interest - coupled to a thermal bath 
responsible for its energy loss. Assuming that any degree 
of freedom of this environment is only weakly perturbed 
by the system of interest, we can describe it as a set of 
independent harmonic oscillators with a coupling which 
is bilinear in the reservoir and system variables and en- 
dowed with a particular spectral function I'l . This model 
has been successfully used to describe general properties, 
classical or quantum mechanical, of dissipative systems 
with only one degree of freedom subject to arbitrary po- 
tentials 1, ^, ^, ^, iQ]. Indeed, it has been extensively 
shown in the literature that, within the range of interest, 
other approaches to dealing with dissipative systems de- 
scribed by a single dynamical variable always furnish us 
with the the same results as those obtained by the bath 
of oscillators with a properly chosen spectral function. 
It is the case , for example, of the application of the 
collective coordinate method to describing the damped 
motion of quantum solitons ^ or microscopic attempts 
to describe more realistic systems such as the electron 
gas of a metallic environment S • 

However, despite all its success there are certain dissi- 
pative systems for which the usual model of the bath of 
oscillators can be shown to be inappropriate to account 
for the physics we expect from them. Here it should be 
stressed that by the usual model we mean independent 
oscillators coupled bilinearly in coordinates to the system 
of interest. 

Suppose one immerses two independent particles in the 
same medium where each of them would separately be- 
have as a Brownian particle. Since for each individual 
particle we could mimic the effect of the medium by the 
bath of oscillators it would be very natural to try to gen- 
eralize the model to cope with the presence of those two 
particles. This generalization is quite straightforward 
and the only point where one should be a bit cautious 



is when introducing the well-known counter-term P, Q 
in the generalized model. In order to do that in an un- 
ambiguous way all one has to do is employ the equivalent 
model for the system-bath Lagrangian where the inter- 
action is described by a velocity-coordinate coupling and 
perform the usual canonical transformations to achieve 
the desired coordinate-coordinate coupling with the ap- 
propriate counter-term 0, I3(see below). It is then a 
simple matter to obtain the equations of motion for each 
particle under the influence of the environment using, 
for example, the Laplace transform method. These are 
coupled Langevin equations that when written in terms 
of the center of mass and relative coordinates, q and u, 
of the two particles present a somewhat bizarre result; 
namely, the latter obeys a free particle equation of mo- 
tion. Therefore, if we give an initial velocity to one of 
those particles the other must instantaneously reply in 
such a way that the relative coordinate evolves linearly 
in time ignoring completely the presence of the environ- 
ment. This is a very unexpected result to say the least. 

It is our intention in this letter to propose an exten- 
sion of the usual model of a bath of oscillators in order 
to fix this deficiency. As we will see, on top of succeeding 
in so doing, we will also be able to describe an effective 
coupling between the two particles mediated by the pres- 
ence of the bath and the resulting interaction potential 
depends on the specific form of the spectral function of 
the environment oscillators. This effect reminds us of 
the formation of Cooper pairs or bipolarons in material 
systems due to the electron-phonon interaction. 

For the sake of completeness we shall start our proce- 
dure by briefly reviewing some general aspects concerning 
the usual model for the bath of oscillators that will be 
useful later on. 

The Lagrangian for the complete system is given by 

Ls + Ln + Li. (1) 

Ls is the Lagrangian of the system of interest, Lr is the 
Lagrangian of the reservoir, which in the case of the bath 
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of oscillators is 
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and Li is the interaction Lagrangian which can be writ- 
ten in two equivalent forms P, Q- For example, in an 
electromagnetic-like fashion we will have a coordinate- 
velocity coupling i/ = — ^ CkRkX. Switching to the 

fc 

Hamiltonian formulation and performing the following 
canonical transformation for the reservoir dynamical 
variables, Pk —>■ mkUJkRk and Rk — > , we can show 
that the coupling term transforms into 



Lr = - 
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(3) 



where the coefficients have been redefined as Ck = Cki^k- 
Therefore, the canonical transformation changes the 
coordinate- velocity coupling into a coordinate-coordinate 
coupling plus a new quadratic term which is necessary to 
preserve the translational invariance of when the sys- 
tem of interest is not acted by any external force. We 
refer the reader to jlOi l for a discussion on the invariant 
form of the latter. 

In our generalization of the previous model, the sys- 
tem of interest with a single degree of freedom will be 
represented by the free particle Lagrangian 



(4) 



The heat bath will be described as a symmetrized col- 
lection of independent harmonic modes, 
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mk(RkR-k - LulRkR- 



(5) 



For the coupling term we initially assume the interac- 
tion Lagrangian 



Lj = -J2Ck{x)Rk. 



(6) 



We have just shown above that this kind of coupling is 
equivalent to the coordinate-coordinate coupling once an 
adequate potential renormalization is introduced. Then, 
following the same procedure that led us to ^ it is easy 
to see the equivalence between (jSJ and 



Li = -IY1 iC-k{x)Rk + Cfc(x)i?_fc) 
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Cfc(x)C-fc(x) 
2mfca;^ 



where Ck{x) = Ck{x)ujk- In order to represent the ef- 
fect of a local interaction of the particle with a spatially 
homogeneous environment we choose 
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With this choice it is easy to show that the entire sys- 
tem is translationally invariant and, moreover, the cou- 
pling to the environment is homogeneous. If the particle 
is displaced by a distance d, the coupling term trans- 
forms into C-k[x + d)Rk = C-k{x)e^'-'''^Rk, which sug- 
gests the definition of a new set of canonical variables as 
Rk — e~'''"^Rk that renders the total Lagrangian invari- 
ant. It is important to notice that with a coupling like 
© , the potential renormalization in Q is a constant and 
therefore does not contribute to the particle dynamics. A 
coupling like lO appears, for example, in the interaction 
of a particle with the density operator of a fermionic bath 
i. 

Now, we write the Euler-Lagrange equations for all 
the variables involved in the problem, take their Laplace 
transforms and eliminate the bath variables, in order to 
find the following equation of motion for the system of 
interest 

Mx+ [ K{x{t)-x(t'),t~t')x{t')dt' ^ F{t). (9) 



Here the nonlinear dissipation kernel, given by 

k "^fe^fc 



K{r, t) = ^ '""'"'2"' cos kr cos cj^r, (10) 



shows that the interaction with the thermal bath induces 
a systematic influence on the system which is non-local 
and non-instantaneous. The function F{t) can be inter- 
preted as a fluctuating force 



C-k{.x)Rk{Q) + Ck{x)R-k{Q) 
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(11) 



where Rk = Rk + '^m'fJs^ ■ ^(^) depends explicitly on the 
initial conditions of the bath variables and its statistical 
properties are obtained from the initial state of the total 
system. Now we need to choose a suitable distribution of 
oscillators, in the continuum limit, which leads us to the 
Brownian motion. The kernel in H1U|) can be written as 
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^7') where lm.x^k\i^) 



cosk{x{t)-x{t'))cosu}{t~t'), (12) 
S{uj — u!k) is the imaginary part 



of the dynamical response Xfc*' {^) of a non-interacting 
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oscillator with wave number k. If we now assume that 
the bath oscillators are in fact only approximately non- 
interacting and replace only their response functions by 
those of damped oscillators, Xk{^), one has 



ImXA;(cj) = 



ruk 



(13) 



where 7fe is the relaxation frequency of that oscillator. 

Our main interest is to study the system for times 
longer than the typical time scale of the reservoir. In 
other words, we wish to study the low- frequency limit of 
eq. H13|l in which Ivaxki'^) oc w and, therefore, we can 
assume 



Imxfc(w) w f{k)uj0{Q, - w). 



(14) 



where f{k) 



Here we have, as usual, introduced 



a high frequency cutoff Q as the characteristic frequency 
of the bath. A functional dependence like l|14() for the dy- 
namical response of the bath has been employed in the 
references for fermionic environments. The partic- 

ular choice of the dynamical susceptibility of the bath al- 
low us to separate its time and length scales and to obtain 
a Markov dynamics when we replace (|14|l in H12|l and in- 
tegrate with respect to lu taking the limit ^ oo. With 
these considerations the equation of motion reads 



Mx{t) +r]x{t) = F{t), 



(15) 



where we have defined rj = k^KkK-kf{k). Notice that 
with this modification we obtain a relation between the 
damping constant and some microscopic parameters of 
the oscillator bath. With the current prescription and 
supposing that the bath is initially in thermal equilib- 
rium it is easy to show that, for high temperatures, the 
fluctuating force F{t) satisfies the relations {F{t)) = 
and {F{t)F{t')) = 2r^kBTS{t - t'), which are character- 
istics of white noise. Notice that this is valid only if we 
assume a classical distribution of oscillators as the initial 
state of the bath. 

In conclusion, the system-reservoir model with non- 
linear coupling presented here allows us to reproduce the 
result one would have obtained by coupling the particle of 
interest bilinearly to a bath of non-interacting harmonic 
oscillators with the spectral function J{lo) — rjuj 1, 2J. 

Now we are going to study the dynamics of a system 
with two degrees of freedom immersed in a dissipative 
environment. In this case the Lagrangian of the system 
of interest is 
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(16) 



and the coupling term 



Li^^IY. ii^-ki^i) + C-k{x2)) Rk 
k 

+ (Ckixi) + Ckix2)) R-k] ■ (17) 



Notice that we have not included any counter-term in 
()17|l since our system is translationally invariant. The 
equations of motion for this Lagrangian are then 

Mx^+ [ K{x,{t)-Xiit'),t-t')i,{t')dt' 



+ / K{x,{t) ~ Xj{t'),t -t')ij{t')dt' 
Jo 

+ -^V{x,{t)-x,it)) = F,{t), (18) 

where i j = 1,2, the fluctuating force Fi{t) has the 
form given by (|ll|l replacing x{t) by Xi(t) and 



V{r{t))^-j:^coskrit). 



(19) 



In the long-time limit, the second term in H18|l can be 
written as r]ii providing us with the usual dissipative 
term. The third term represents a cross-dissipative term 
that depends on the velocity of the second particle and 
the relative distance between them. The last term is 
clearly an effective interaction induced by the coupling 
with the thermal reservoir. To see the explicit form of 
this interaction and the non-local influence in the cross- 
dissipative term, we need to evaluate the sum in H12() 
with lTaxf\i<j) replaced by (|14|l . For this we had better 
define the function g(k) as 



Zn 



(20) 



where we have assumed a one-dimensional environment 
of length L. 

The exact form of this function can be obtained from 
a microscopic description of the bath and the details of 
the interaction between the system of interest and the 
bath degrees of freedom. As we are not interested in 
such a detailed description of the many-body reservoir, 
we are going to choose a form for g{k) that satisfies the 
condition g{k)k'^dk — 1 derived from the definition 
of the damping constant r). Notice that we have replaced 
^ J dfc. A reasonable choice for this function is 



(21) 



where A = 1/(2A;q) is the normalization constant and fcg 
determines the characteristic length scale of the environ- 
ment. For example, in the case of fermionic environments 
fco is of the order oi kp |^. One should bear in mind that 
the choice of g{k) is guided either by the knowledge of 
the microscopic details of the environment or by some 
phenomenological input about the effective interaction 
between the two particles. 

With the form chosen for g{k) the kernel of the third 
term in H18|l is, again only for high temperatures. 



K{xi{t) - X2{t'), t-t')^ 2S{t - t')r]{u{t)), (22) 
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Figure 1: Spatial Correlation of -Fi and F2 



where u{t) = xi{t) — X2{t) and 

\ (/cjfii^ _|_ 1) (/cju^ + 1) y 

and the effective potential reads 

The strength of the effective potential depends on the 
characteristic length and time scales. Therefore the con- 
tribution of this coupling to the dynamics of the Brown- 
ian particles is only important for times longer than the 
typical time scale of the reservoir and distances of the 
order of (or less than) the characteristic length k^^. The 
fluctuating forces still satisfy the white noise properties 
but now they present an additional property associated 
with the distance between the particles. The forces Fi{t) 
and F2{t) are also spatially correlated, that is 

{Fiit)F2it')) ^2Tj{u{t))kBTS{t - t'). (25) 

In figure ^ we see the noise correlation strength as a 
function of the distance between the particles. For short 
distances the noise correlation has the standard Brown- 
ian behavior. However, for longer distances the correla- 
tion function becomes negative and this anti- correlation 
induces an anomalous diffusive process in the system 
which will ultimately tend to normal diffusion once the 
particles are infinitely far apart. 

In terms of the relative and center-of-mass coordinates 
the equations of motions read 

Mu{t) + r]u{t) - 7j{u{t))u + V'{u{t)) = Fu{t) (26) 

and 

Mq{t)+m{t)+ri{u{m = Fq{t), (27) 

where i^„(t) = Fi{t)-F2{t) &ndFq{t) = (Fi (i)+F2(i))/2. 
^From the form of V{u{t)) and rj{u{t)) and the sta- 
tistical properties of the fluctuating forces, {Fu{t)) = 



{Fq{t)) = 0, {Fu{t)Fu{t')) = AkT{rj - v{u))S{t - t') and 
{Fq{t)Fq{t')) = kT{ri + 'n{u))S{t - t'), it is evident that 
at large distances the equations of motion for the rela- 
tive and center of mass coordinates represent the motion 
of Brownian particles with reduced mass M/2 and total 
mass 2M, respectively. In general the dissipation de- 
pends on the relative distance between the particles and 
for distances such that uko <C 1, we have up to first order 
in u(t), V'{u) oc —u{t). In this approximation both dis- 
sipation and fluctuations are negligible and then we have 
an undamped oscillatory motion for u{t). 

In conclusion, we presented a system-plus-reservoir 
model with a coupling which is non-linear in the system 
coordinates and, in the adequate limit, allows us to repro- 
duce the phenomenological results known for the dynam- 
ics of a dissipative system with only one degree of free- 
dom. In our model the dissipation coefficient is directly 
expressed in terms of a few parameters of the thermal 
bath that can be measured experimentally. Moreover, 
our model is capable of inducing an effective coupling 
between the two Brownian particles which arises from 
the non-local effects generated by the bath and depends 
on the choice of the dynamical susceptibility of the reser- 
voir. In the single particle case, the non-local effects will 
be relevant only in a non-markovian approximation. 

We should finally stress that we have only reported 
our findings for the classical or high temperature regime. 
The quantum mechanical analysis of this system requires 
the evaluation of the reduced density operator for the two 
particles of interest down to very low temperatures and 
will be presented elsewhere [lH . 
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